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I consider the problem of vortex tunneling in a two-dimensional superconductor. The vortex 
dynamics is governed by the Magnus force and the Ohmic friction force. Under-barrier motion in 
the vicinity of the saddle point of the pinning potential leads to a model with quadratic Hamiltonian 
which can be analytically diagonalized. I find the dependence of the tunneling probability on the 
normal state quasiparticle relaxation time r with a minimum at loqt ~ 1, where u>o is the level 
spacing of the quasiparticle bound states inside the vortex core. The results agree qualitatively with 
the available experimental data. 

PACS numbers: 74.60. Ge 



I. INTRODUCTION 

The dynamics of vortices in Type-II superconductors 
has been one of the most attractive areas of research in 
recent years™ because of its importance from scientific as 
well as technological point of view. To science, vortices 
present an example of extended topologically stable ob- 
jects with extremely rich properties. On the technological 
side, the motion of vortices is the source of supercurrent 
dissipation - the circumstance that limits high magnetic 
fields application of high-T c materials. In view of this 
fact, the exploration of the limits of strong pinning of 
vortices has been actively pursued both in theory and 
experiment. It has long been known that the thermally- 
activated depinning rate strongly decreases as the tem- 
perature is lowered and eventually saturates at values 
believed to be set by quantum tunneling. The present 
article explores this quantum regime. 

The difficulty of the problem is in the fact that there 
are two major forces that govern the dynamics of the 
vortex: Hall (Magnus) force and friction force. While 
the first one is conservative and can be treated through 
a single vortex description, the second one is not: the 
energy dissipates into the environment. Therefore, one 
needs to consider the combined vortex-plus-environmcn± 
system. In the pioneering paper by Feigel'man et al.fi 
which took into account both forces, the environmental 
degrees of freedom were integrated out with the use of 
path-integral techniques to produce an effective vortex 
action description. The time non-locality of the resulting 
action, however, limited the analysis to qualitative con- 
clusions. Later, the same approach was independently 
taken by Morais Smith et alp but they also failed to go 
beyond the scaling analysis of the effective action. 

This paper is devoted to the study of an exactly solv- 
able model of two-dimensional dissipative vortex tunnel- 
ing. It has a quadratic Hamiltonian of the vortex coupled 
to the environment and is solvable by an analytic diago- 
nalization. The present model is a generalization of the 
quadratic model of one-dimensional dissipative particle 



tunneling that was solved by Ford et a/B The results of 
the solution show the dependence of the tunneling rate 
on the ratio of the Magnus and dissipative forces. The 
tunneling rate has a minimum near the point where these 
two forces become equal. 

This paper is organized as follows: I first start with 
a non-dissipative case in Sec. 2; the dissipative model is 
introduced in Sec. 3, where some of its general features 
are analyzed; in Sec. 4 I describe the solution of the dis- 
sipative tunneling problem and compare the results with 
available experimental data; conclusions follow in Sec. 5. 
Units with Ti — 1 are used throughout the paper. 



II. VORTEX TUNNELING: NON-DISSIPATIVE 

CASE 

In order to analyze the motion of a vortex one needs 
first to identify the forces that govern its dynamics. Both 
in classicalO and in quantumQ fluids there is an intrinsic 
Magnus force that acts on a vortex when it moves rela- 
tive to a (super)fluid. This force is normal both to the 
vortex tangent and to the vortex velocity, relative to the 
fluid, and is linear with in the latter. In this respect, it is 
very similar to Lorentz force that acts on charged parti- 
cles moving in magnetic fields. Since Magnus force does 
not produce any work, it is possible to fownulate vortex 
dynamics using Hamiltonian formalism]!™ This formal- 
ism was successfully applied to the problem of quantum- 
mechanical vortex nucleation by Volovik.Q Howeyer, it 
was not until much later that Haldane and Wvj3 real- 
ized that, just like Lorentz force is ra manifestation of 
a geometric (Aharonov-Bohm) phaseS associated with a 
motion of a particle, Magnus force should arise as a con- 
sequence of a geometric phase-associated with vortex mo- 
tion. Later, Ao and ThoulessE3 extended this idea to vor- 
tices in superconductors. To clarify the physics involved, 
I will sketch the derivation of the geometric phase done 
by Ao and Thouless. 
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Consider the ground state many-body wave function 
of a two-dimensional superconductor: fo( r ivi r w), 
where rj = (xj,yj) are the coordinates of electrons. A 
vortex at r = (x,y) can be created by "spinning up" the 
system: 



. N 

.7=1 



(r 3 -r) } * (ri,...,rjv;r). 



(1) 



Here it is assumed that ^ i s "close" to *■ e., it is 
obtained from the latter by a continuous deformation; 

6{rj — r) = arctan ( ^/.Zy ^J ■ Note that the factor i in 
the exponent comes from the fact that the condensate is 
made of Cooper pairs. Now one can adiabatically move r 
around a closed contour keeping the coordinates of elec- 
trons fixed. The phase accumulated at the end will come 
from the exponent and will be given by tt times the num- 
ber of electrons encircled by the path of the vortex. This 
phase can be described by the inclusion of a geometric 
phase term in the action: 



/ 



dtax(t)y{t). 



(2) 



Here a = irn s is the Hall coefficient and n s is the density 
of electrons in the condensate. One sees that the length 
scale, associated with this geometric phase term is 1/ ^/n^ 
which is one of the smallest in the problem. One expects 
then that the dynamics will be dominated by the Magnus 
force. 

Next I include the vortex potential energy in the ac- 
tion: 



S = dt [axy - V(x,y)] 



(3) 



To clarify the picture I go over to the Hamiltonian for- 
mulation. First of all, it is easy to identify canonically 
conjugate variables. Indeed, from the above action one 
sees that the variable, conjugate to x is ay. The Hamil- 
tonian is then equal to the potential energy: 



[x,y] 



H = V(x,y), 
i 

a 



(4) 
(5) 



For an arbitrary V{x,y) there is a problem of operator 
ordering. i-This problem has been addressed by Girvin 
and Jach;Eil an alternative approach to this problem is 
the path-integral approach considered by Jain and Kivel- 
son.E3 In the model studied below, however, this difficulty 
does not appear. 

A comment should be made about the relevance of vor- 
tex mass which is not included in the model. Using the 
analogy between a superconducting vortex and a particle 
in magnetic field (i. e. similar origins of the correspond- 
ing geometric phases) this omission of the vortex mass 




FIG. 1. Saddle point S in the vortex pinning potential bi- 
ased by a supercurrent J (arbitrary units). 



amounts to the effective Lowest Landau Level approxi- 
mation. Such an approximation is warranted if the tran- 
sitions to higher Landau levels can be neglected. This can 
be checked by comparing matrix elements |V^ m | of the 
potential energy V(x, y) between different Landau levels 
with the spacing oj c between Landau levels. The ratio 
I V nm | /uj c depends on the particular potential under con- 
sideration and can be estimated in our problem of vortex 
tunneling to be ~ 2-k / S, where S is the action involved 
in tunneling. The value of this action can be extracted 
from experimental data (to be discussed in Sec. 4) and is 
w 50. Thus, the ratio \ V nm \/ui c ~ 0.1 and the omission of 
the vortex mass is indeed warranted. However, it is pos- 
sible that in some other situations the mass of a vortex is 
relevant and should be included. This can be easily done 
since the mass term in the Hamiltonian is quadratic and 
thus does not lead to any substantial complications. 

Consider now the tunneling problem in which an im- 
purity potential, which pins the vortex, is tilted by a 
Magnus force due to an applied depinning supercurrent 
(see Fig. |l|). For a vortex to tunnel out of the pinning site 
it has to overcome a potential barrier. This leads to a 
thermally- activated depinning rate at high temperatures. 
Here I am considering the problem of quantum tunnel- 
ing which gives dominant contribution to the depinning 
rate at low temperatures. One can argue then that for 
sufficiently strong depinning supercurrents the tunneling 
exponent is dominated by the under-barrier motion in the 
vicinity of the saddle point of the potential (see Fig. [l]). 
This leads one to consider the problem of tunneling in 
the following model potential: 



= Ky 2 _ 2 

2 y 2 ' 



> 0. 



(6) 



This becomes the problem of tunneling across an inverted 
parabolic barrier — k x x 2 in ID quantum mechanics if one 
identifies the momentum p — ay, and the mass m = 
a 2 1 ' n y of the ID particle. The transmission coefficient 
forJJiis problem can be calculated exactly and is given 
byB 
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D (E) = 



1 



1 + cxp 



(7) 
(8) 



Here — E is the energy of the vortex as measured from 
the value at the saddle point. E equals the activation 
energy in the high temperature thermally-assisted tun- 
neling regime. 

In comparison, the previous works on the dissipative 
vortex tunneling problem&tl considered the model with 
an added cubic term: V = ^-y 2 — ^x 2 — ^x 3 which 
has a minimum at x = — — . I argue that the role of the 
cubic term is mainly to set the position of the localized 
state and its energy and that it is not essential for the 
tunneling itself. Meanwhile, setting K3 — produces a 
quadratic Hamiltonian which allows for an exact diag- 
onalization when the dissipation is added. Finally, let 
me mention that the model Eq. (|^) (without dissipation) 
has been considered before by Fertig and HalperinlU in 
a different context: the tunneling of electrons in two di- 
mensions in the presence of a strong magnetic field. It 
later became the basis for the Chalker-Cpddington model 
of the Integer Quantum Hall transition.ll3 

The result Eqs. (0,||) shows an interesting feature of 
the vortex tunneling problem. If n y is decreased while 
both k x and the height of the barrier E are kept con- 
stant the transmission coefficient Dq decreases implying 
stronger pinning. In particular, for n y = the vortex no 
longer tunnels across the barrier. This result can be un- 
derstood by observing that for K y = the Hamiltonian 
becomes invariant with respect to the y coordinate of the 
vortex. This leads to the conservation of the conjugate 
momentum which is ax. Thus, the tunneling process 
in which x changes is forbidden by a new conservation 
law. This means that the pinning potentials which are 
translationally invariant along the direction of the bias- 
ing supercurrent are good candidates for strong vortex 
pinning. 



III. DISSIPATION 

So far we have considered a vortex in a pure supercon- 
ductor. In reality, the presence of impurities that sets the 
normal state resistance at low temperatures introduces a 
finite quasiparticle scattering time r. This fact has to 
be taken into account and the adiabatic arguments pre- 
sented above need to be modified. These modifications 
have been worked out b#, Kopnin and KravstsovE^I and 
by Kopnin and SalomaacH (for a simple new look at the 
problem see alsotj). They showed that a finite relax- 
ation time breaks the adiabaticity in the spectral flow of 
the quasiparticle vortex core bound statealj and leads to 
two main effects: the Hall coefficient a is reduced from 
its pure value im s , and there appears a friction force act- 



ing on the vortex, F 
given by: 

M 2 



-r/v. The coefficients a and r\ are 



' 1 + (u; r) 2 



n 



UJ T 



1 + (loot) 2 



(9) 



Here, ojq is the level spacing of quasiparticle bound states 
inside the vortex core. Thus, the parameter ujqT sets rel- 
ative importance of the Magnus and friction forces. This 
can be quantified by defining the Hall angle: tan0# = 
a/r/ = loot. Thus, only in the "superclean limit" Uqt ^> 1 
can one neglect the dissipative force and consider a pure 
quantum problem. 

It should be mentioned that these is an alternative 
point of view expressed by Ao et alS-zi according to which 
the Hall coefficient is a topological number and thus is 
not renormalized. It turns out that-the analysis of exper- 
imental results of van Dalen et alSB presented in Sec. 4 
suggests that there is indeed a renormalization of the 
Hall coefficient a. Meanwhile, I take a and 77 as two 
phcnomenological parameters. n 

In the following I use the approach of Ford et ala to the 
problem of dissipation in quantum mechanics. I model 
the friction force by a linear coupling to a bath of oscil- 
lators (to simplify notation I use x a , a = {1, 2} for x and 
y coordinates of the vortex) : 



H = V(x a ) + J2 



2 m 1 



1 



(10) 



[X1,X 2 ] = -, [qai,Pbj] = iSabStj. (11) 

a 

Vortex coordinates x\ and xi are coupled to two inde- 
pendent identical sets of oscillators in a simple way: they 
shift their equilibrium positions. The difference between 
this model the one considered by Ford et ala is that, due 
to a different physical context, here both x\ and x-i need 
to be coupled to the oscillator baths in order to account 
for friction in both directions. 

The Heisenberg equations of motion corresponding to 
Eq. @ are: 



ae ab Xb 



= ~ dx~ + ^ mj j ^ aj ~ Xa '' 
j 

q aj + Lo]q a] = uj 2 x a . 



(12) 
(13) 



The solution of Eq. (13) is: 



Qaj(t) = q h a At) + x a (t) - / alt' cos[^-(t - t')]x a (t'), (14) 



where q^ (t) is the solution of the corresponding homoge- 
neous equation. Substituting Eq. ( |l4| ) into Eq. ( |l2| ) one 
obtains: 
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dV 

ae ab x b = -— - I dt' fi(t-t')± a (t')+N a (t), (15) 



//(i) = 0(f) m j ljj2 j cos(u>jt), 

3 

N a (t) =^^(4). 



(16) 



(17) 



Here //(f) is the so-called memory kernel, N a (t) is the 
noise force and 0(f) is the Heaviside step function. 

Next I consider the so-called Ohmic case when the fric- 
tion force acting on the vortex is strictly linear in the 
vortex velocity. This is realized by taking the following 
distribution of oscillator strengths: 

rjioS) = — m j UJ 'j [S((*> — &j) + S(uj + LUj)] = const. (18) 
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FIG. 2. Magnetization relaxation rate Q (normalized by its 
value Qoo at luqt — > oo) vs. loot plotted for various values of 

S = K x /Hy ~ yjl — J/Jc- 



With this choice Eq. (|15|) becomes: 



9x a 



??.T a +7V a (f). 



(19) 



Naturally, this can be called the force balance equation: 
it balances Magnus, potential, friction and noise forces. 
This equation shows that the model Eq. ( |To| ) indeed de- 
scribes friction force linear in the velocity. A more gen- 
eral version of this., equation is known as the quantum 
Langevin equation.!] 



IV. DISSIPATIVE VORTEX TUNNELING 

Returning to our problem of tunneling in the vortex 
depinning by an external supercurrent, one now has to 
consider Eq. ( [To|) with V given by Eq. (||). The cru- 
cial observation due to Ford et aln (made for a similar 
problem) is that the resulting system is a set of coupled 
oscillators with one of them having a negative spring con- 
stant (— k x ). This leads to a spectrum which contains an 
isolated mode with purely imaginary frequency ifl*. Af- 
ter the diagonalization of the Hamiltonian one obtains an 
oscillator with a (shifted) negative spring constant that 
is decoupled from the rest of the oscillators (all with pos- 
itive spring constants). It is this new negative spring 
constant that determines the tunneling amplitude. The 
transmission coefficient is then given by Eq. (Q) with a 
substitution f^o — ^ 



D(E) 



1 



l + exp{2^} 



(20) 



The value of the imaginary eigenmode frequency iCl* can 
be found by diagonalizing the Hamiltonian; a much eas- 
ier way is to observe that this ei gen mode should satisfy 
the force balance equation, Eq. ([19]), without the noise 
term (i. e., the corresponding equation for the expec- 
tation values; it should be stressed that noise term can 



be separated due to the linearity of the equations) . This 
leads to the following equation on O*: 



det 



— k x + r]Q* 



Ky + ryfT 



= 0. 



(21) 



It is easy to understand this equation qualitatively us- 
ing the following simple arguments: roughly, all oscilla- 
tors in the bath (see Eq. (|l0|)) can be divided into two 
types, those with frequencies smaller and bigger than Q* . 
Oscillators with small frequencies do not follow the mo- 
tion of the vortex and each one effectively increases the 
values of n y and —k x by rrijUXj. Those with high fre- 
quencies adjust to the vortex motion and thus do not 
contribute to the increase of k's. Then a simple estimate 
gives : 



n y 

- h :r 



K y + £ rrijw] 



<<n* 



- K x 



Crin*, 



with C ~ 1. This indeed agrees qualitatively with the 
correct Eq Q . 

Eqs. ( p0||2l| ) constitute the solution of the problem. In 
the rest of this section the implications of this solutioa 
for the existing experimental data of van Dalen et a/.CJ 
will be discussed. In these experiments quantum dynam- 
ical relaxation rate of magnetization Q was measured in 
oxygen depleted thin films of YBa 2 Cu30 x as a function 
of the oxygen content x. Changing x changes the normal 
state quasiparticle relaxation time r which, in turn, is 
expected to change in the Hall and friction coefficients. 
I use expressions in Eq. (||) to compare the results of the 
solution of the present model with the experiment. 

The dynamical zero temperature magnetization relax- 
ation rate is Q = 1/5 where 5 is the action involved in 
a tunneling event. It determines the transmission coeffi- 
cient of the pinning barrier as D ~ exp(— 5). In reality Q 
is small: Qoo = Q{ujqt — ► oo) ps 0.02, therefore Eq. ( p0| ) 
implies: Q — Q,* /2nE. As the oxygen content, and thus 
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ujqt, is varied fT (given by the solution of Eqs. (pl|Js|)) 
changes while E (the height of the pinning barrier) stays 
constant. In Fig. || the calculated ratio Q/Qoo is plotted 
as a function of ujqt for several values of 6 — k x /k v . A 
simple calculation gives that S ~ ^/l — J/J c , where J c is 
the critical biasing current density at which the minimum 
of the pinning potential disappears. 

From Eq. (g) one can see that the Hall angle is 
tan Qh — ct/rj = luqt, therefore in the "superclean limit" 
luqt — ► oo the dynamics is dominated by the Magnus 
force, while in the opposite "dirty limit" ujqt — > it is 
dominated by the force of friction. Fig. |^ shows that in 
the superclean limit Q saturates while in the dirty limit it 
varies as 1/{u>qt). Qualitatively, both of these behaviors 
were known before.Em What was not known, however, 
was the behavior of Q in the crossover region luqt ~ 1. 
The new qualitative feature of the exact solution of the 
present model is a minimum of Q in the crossover region; 
it is more pronounced for stronger depinning currents. 
This result implies that the strongest pinning occurs in 
the "moderately clean" regime when ujqt is pL order 1. 

In the experiments of van Dalen et alEi only the 
regime ujqt < 1 was realized. There Q shows roughly 
a 1/(loqt) dependence. In the region luqt ~ 1 the onset 
of a crossover in Q is visible, but the minimum of Q can 
hardly be observed since no data are available in the su- 
perclean limit ujqt 3> 1. In view of this, experiments on 
cleaner films would be desirable to check the existence of 
the minimum in Q - the main qualitative predictions of 
the above analysis. 

Finally, it should be emphasized that the 1/{ujqt) be- 
havior of Q in the dirty limit that was observed in the 
experiment is a clear evidence in favor of the correct- 
ness of Eq. (^J). Indeed, if the Hall coefficient a were not 
renormalized by disorder it would imply that Q in the 
dirty limit should always be smaller than its superclean 
limit value Qoo due to ±he general tendency of dissipation 
to suppress .tunneling. E3 The experimental results of van 
Dalen et alxB show that this is not the case. 



V. CONCLUSION 

In this article, I have considered a model of dissipa- 
tive vortex tunneling. In this model the system-plus- 
environmcnt Hamiltonian is quadratic and can be ana- 
lytically diagonalized. The results were obtained for the 
dynamical magnetization relaxation rate Q as a function 
of the Hall angle tanO# = a/r/ = luqt, where luq is 
the level spacing of the quasiparticle vortex core bound 
states. The results show a 1/(u) t) dependence of Q in 
the "dirty" limit ojqt <C 1, saturation in the "superclean" 
limit ujqt 1 and a minimum at lo^t ~ 1. These pre- 
dictions were compared .with the available experimental 
data of van Dalen et a/..E3 Results agree qualitatively in 
the "dirty" regime ujqt <C 1. On the other hand, the 
predicted minimum can not be found since no data are 



available in the "superclean" regime ujqt ^> 1. 

Although I have only considered here the case of a vor- 
tex in a two-dimensional film, the results can be carried 
over to the three-dimensional case as well. Following the 
arguments of Brandt ,E3 one can show that the effect of 
the extra elastic term in the energy of a vortex line is 
mainly to set a length scale along the vortex - the size 
of the tunneling nucleus. Meanwhile, the effect of dis- 
sipation on vortex tunneling should be qualitatively the 
same. 
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